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In X-ray interference fringes accompanied by mirage diffraction, variations have 
been observed in the spacing and position of the fringes from a plane-parallel Si 
single crystal fixed at one end as a function of distance from the incident plane of 
the X-rays to the free crystal end. The variations can be explained by distortion 
of the sample crystal due to gravity. From the variations and positions of the 
fringes, the strain gradient of the crystal has been determined. The distribution 
of the observed strain agrees with that expected from rod theory except for 
residual strain. When the distortion is large, the observed strain distribution 
does not agree with that expected from rod theory. 



1 . Introduction 

Under Bragg geometry, a refracted beam of X-rays in a weakly 
bent crystal propagates along a hyperbolic trajectory and 
returns to the incident surface, and then the diffracted beam is 
emitted from the surface (Gronkowski & Malgrange, 1984), as 
depicted in Fig. 1. The beam corresponding to the Poynting 
vector of the X-ray is called the refracted beam in this paper. 
Since Authier (2001) has pointed out that the behavior of the 
refracted beam in a distorted crystal resembles a mirage of 
light, we call this diffraction beam a mirage diffraction beam. 
The Interference Fringes between the Mirage Diffraction 
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Figure 1 

Beam trajectories in a weakly bent crystal and schematic illustrations of 
interference between the beams, (a) Interference of a mirage diffraction 
beam (dashed curve with the corresponding refracted beam Si) with 
another mirage diffraction beam (solid curve with the corresponding 
refracted beam S 2 ) to produce IFMD in a weakly bent part of the crystal. 
(b) Interference of a mirage diffraction beam (solid curve with the 
corresponding refracted beam S 3 ) with a reflected beam from the bottom 
surface (dashed curves with the corresponding refracted beam S 4 ) to 
produce IFMRB in a more weakly bent part, h is the reciprocal lattice 
vector and u the atomic displacement vector, and a l , ot 2 , ot 3 and a 4 are the 
glancing angles of the incident beam corresponding to the refracted 
beams Si, S 2 , S 3 and S 4 , respectively. 



beams Si and S 2 in Fig. 1(a) (IFMD) may be observed in the 
beam emitted from A 2 . IFMD were first observed from a 
quartz crystal by Zaumseil (1978). More recently, IFMD were 
observed from an Si strip crystal that was fixed at one end and 
distorted by applying a force vertically at the other end (Yan et 
al, 2007; Fukamachi et al, 2010; Jongsukswat et al, 2012a). 
The IFMD fringe spacing becomes small as the distance (x) 
between the incident plane of the X-rays and the emission 
point of the diffracted beam increases, or as the strain gradient 
increases. The strain gradient s' zz of the strip was measured by 
changing the bending displacement (D) at the point of action 
on the strip. The determined values of s' zz were in good 
agreement with those calculated from rod theory (Jong- 
sukswat et al, 2012a). When the crystal distortion becomes 
very small and the IFMD disappear, another set of Inter- 
ference Fringes between a Mirage diffraction beam S 3 and a 
beam Reflected from the Bottom surface S 4 in Fig. 1(b) 
(IFMRB) was observed (Fukamachi, Jongsukswat et al, 2011). 
The IFMRB fringe spacing becomes large as the distance x 
increases, which is the opposite variation to the IFMD spacing. 

In this paper, we report that we have observed both IFMRB 
and IFMD interference fringes from a plane -parallel Si single 
crystal fixed at one end when no artificial force is applied at 
the free end. The force of gravity applies uniformly to the 
sample in the present case, which is different from the force in 
the previous experiments (Fukamachi et al, 2010; Jongsukswat 
et al, 2012a). We have compared the measured values of s' zz 
with those calculated from rod theory by assuming a uniformly 
distributed load. 

2. Experimental 

In the experiment, two plane-parallel Si single crystals were 
used. The two sample crystals were cut from the same wafer 
and prepared in the same way. The top and bottom surfaces of 
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Figure 2 

(a) A photograph of the sample fixed to the holder, (b) Sample and X-ray beam geometries. L is the 
distance from the fixed end to the free end, / the distance from the free end to the incident plane of 
the X-rays and D the displacement of the crystal end from the flat surface, (c) A schematic diagram 
of the measuring system. x L is the distance between the incident point of the X-rays and the crystal 
edge in the incident plane. 



the crystals were polished at Sharan Instruments Corporation 
using a non-disturbance polishing technique. Both samples 
were 15 mm wide and 0.28 mm thick. One end of the first 
sample was pasted to an Al plate with pine resin, as shown in 



Figs. 2(a) and 2(b). The distance 
between the free and fixed ends of the 
sample was 41 mm. The second sample 
was used in the measurement described 
in the latter part of §3. 

The beam geometry around the 
sample is shown in Fig. 2(b). It is noted 
that the incident plane of the X-rays is 
normal to the bending direction (the 
longer dimension) of the sample. A 
schematic diagram of the optical 
measuring system is shown in Fig. 2(c). 
Si(220) diffraction experiments were 
carried out using X-rays from 
synchrotron radiation at the bending 
magnet beamline BL-15C, Photon 
Factory, KEK, Tsukuba, Japan. The 
X-rays were er-polarized and mono- 
chromated by an Si(lll) double-crystal 
monochromator. The X-ray energy was 
11 100 eV, with the error in the spectral 
center being 0.5 eV. The vertical and 
horizontal widths of Slit 1 were 0.02 
and 4 mm, respectively. In Fig. 2(c), P h 
represents the intensity of the 
diffracted beam (primary peak) from 
the incident point, P 1 ^ represents that 
of the mirage diffraction beam from 
the surface and P^' that of the emitted 
beam from the lateral surface in the 
direction of the diffracted beam. The superscript (n) repre- 
sents the serial number of the fringes. The beam intensities 
were recorded on a nuclear plate (Ilford L4, emulsion thick- 
ness 25 um) and measured using a scintillation counter (SC). 
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Figure 3 

A series of section topographies from an Si plane-parallel crystal in the Bragg mode, taken by changing the position of the incident plane (horizontal 
direction corresponding to /). The vertical direction from the top corresponds to the direction of x. x L = 7.5 mm in the right-most panel. 
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Fig. 3 shows a series of reflection topographies of Si(220) 
recorded by changing the distance / from the free end to the 
incident plane of the X-rays (from right to left in Fig. 3) in 
Fig. 2(b). The incident plane of the X-rays was varied from the 
free end to approximately the fixed end. The dark contrast in 
the uppermost part of Fig. 3 is the primary peak of P h , while 
the darker contrasts in the lowest part are the interference 
fringes observed in the emitted beam from the lateral surface 
in multiple Bragg-Laue mode, P^ (Fukamachi, Hirano et al., 
2011). The interference fringes observed between P h and P^' 
are either IFMD or IFMRB. Since the spacing of the inter- 
ference fringes increases as x (the distance from the top in the 
vertical direction in Fig. 3) increases for / < 12 mm, these 
interference fringes are attributed to IFMRB. On the other 
hand, since the spacing of the fringes decreases as x increases 
for / > 24 mm, the interference fringes are attributed to IFMD. 
In the intermediate region 12 < I < 24 mm, both IFMD and 
IFMRB are observed. The spacing of the interference fringes, 
either IFMD or IFMRB, decreases as the strain or / increases, 
as can be seen in the variations along the horizontal direction. 

3. Analysis of interference fringes 

In order to analyze the variations in spacing of IFMD and 
IFMRB, we use the dynamic theory of diffraction for a weakly 
distorted crystal, because we know that such variations are 
caused by distortion in the crystal (Fukamachi et al, 2010). On 
the basis of the dynamic theory of X-ray diffraction in Laue 
mode for a weakly distorted crystal developed by Penning & 
Polder (1961), Gronkowski & Malgrange (1984) extended the 
theory to the symmetric Bragg mode, except for the total 
reflection region. According to the theory, the beam trajectory 
of the refracted beam of mirage diffraction is given by 



+ W s j - [fix + s(W s )(Wt - l) 1/2 f = 1, (1) 

where the direction of x is taken parallel to the crystal surface 
and that of z in the direction of crystal depth, as shown in 
Fig. 1(b). The value of s(W s ) is 1 for W s > 0 and -1 for W s > 0, 
where W s is the value of the deviation parameter W of the 
incident X-ray at the incident point. The parameter W is 
defined by 
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(2) 



Here, a is the glancing angle, 6 B the Bragg angle, C the 
polarization factor, and Xo and Xh the 0 and hth Fourier 
components of X-ray polarizability, respectively. According to 
Yan et al. (2007), the parameter fi is given by 



P = 



2 sin 2 6 B tan 6 n 



C(XhX-h) 



1/2 £ zz ' 



(3) 



where e' zz is the strain gradient. Under an anomalous trans- 
mission condition (W s < —1) and with Wfi < 0, the beam 
trajectory shows a hyperbolic form, as shown in Fig. 1, and the 
eccentricity is related to /3. 



According to Yan and co-workers (Yan & Noyan, 2006; Yan 
et al, 2007), the strain e zz of a plane -parallel crystal is given by 



=^='-('-D- 



(4) 



Here, d is the spacing of the diffraction planes under stress and 
d 0 that without stress, and H is the crystal thickness. The value 
of z is taken to be 0 at the top surface and H at the bottom 
surface. According to rod theory, e' zz may be written as 

4vD 



<z=Sil + A<4 withg ; = 



L 4 



(5) 



when the crystal is fixed at one end and distorted by a 
uniformly distributed load due to gravity. In equation (5), v is 
the Poisson ratio and As' zz the residual strain gradient. In an 
ideally plane-parallel crystal, As' zz should be zero. 

Because IFMRB are observed for / < 14 mm in Fig. 3, a 
value of s' zz in this region is obtained by fitting the IFMRB 
profile calculated by assuming a value of j3(e' zz ) to the 
measured IFMRB profile for each value of I, as described by 
Fukamachi, Jongsukswat et al. (2011) and Jongsukswat et al. 
(20126). The insets of Fig. 4 show an example of the fit. The 
measured IFMD and IFMRB profile for I = 14 mm is shown in 
the upper inset, and the best-fit curve calculated by taking ten 
mirage diffraction beams into account is shown in the lower 
inset. On the other hand, IFMD are observed for / > 14 mm. 
The value of s' zz in this region is obtained by measuring the 
position of the third peak for each value of /, as described by 
Jongsukswat et al. (2012a). It is confirmed that the value of s' zz 
is also obtained by fitting the position of the third peak to the 
measured one, and this value is the same as that obtained by 
measuring the peak position. The measured position of the 
third IFMD peak is indicated by an arrow in the upper profile 
and the calculated one is indicated by an arrow in the lower 
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Figure 4 

The measured values of e' zz as a function of I 2 . The solid circles show the 
values obtained using IFMRB for I 2 < 196 mm 2 and IFMD for I 2 > 
196 mm 2 . The line shows the best fit to the measured values. The insets 
show an example of the best-fit calculated IFMD and IFMRB curve to the 
measured profile for / = 14 mm. The position of the third peak of the 
measured IFMD is indicated by an arrow in the upper profile and the 
corresponding calculated position is indicated by an arrow in the lower 
curve. 
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Figure 5 

(a) A photograph of the cantilever for applying force in the opposite direction to gravity, (b) The geometry of the sample and X-rays, (c) A schematic 
diagram of the measuring system after the Si(lll) double-crystal monochromator in Fig. 2(c). An Si(220) perfect crystal is inserted between the double- 
crystal monochromator and Slit 1 in Fig. 2(c). Slit 1 is rotated 34° anticlockwise and Slit 2 34° clockwise from the positions in Fig. 2(c). The vertical widths 
of Slit 1 and Slit 2 are both 40 um. 



curve. The solid circles in Fig. 4 show the resultant values of 
the strain gradient (s' zz ) as a function of I 2 , obtained using 
IFMD or IFMRB in Fig. 3. The values of e' zz in Fig. 4 show a 
linear increase as a function of I 2 in the range 0 < I 2 < 800 mm 2 , 
which agrees with rod theory. The fitted line is also shown in 
Fig. 4. The value of g t in equation (5) is obtained from the 
slope of the line and the value of As' zz is obtained from the 
intercept of the line with the line of I 2 = 0 mm 2 . The values 
thus obtained are g t = 4.9 x 10~ 8 mm~ 3 and As' zz = 3.8 x 
10~ 6 mm -1 . By putting the obtained value of gt = 4.9 x 
10~ 8 mirT 3 , together with v = 0.279 and L 2 = I 2 ~ 800 mm 2 , 
into equation (5), the displacement D is determined to be 
28 um. In the range I 2 > 800 mm 2 , the observed values of e' zz 
become larger than those expected from the linear depen- 
dence. This discrepancy may be explained by taking the 
pasting force to the sample holder into account in rod theory. 

In the experiment in Fig. 3, no force was applied artificially 
at the free end of the sample and the strain was caused by the 
force of gravity uniformly loaded on the sample. In order to 
confirm this, we applied a force in the opposite direction to 
gravity. A photograph of the cantilever is shown in Fig. 5(a), 
and the geometry of the beam arrangement and the applied 
force is shown in Fig. 5(b). The measuring system after the 
Si(lll) double-crystal monochromator in Fig. 2(c) is shown in 
Fig. 5(c). The second sample was used. The length L of the 
sample was 40 mm. The sample was held between two sheets 
of aluminium fixed with bolts at both ends of the sheets, as 
shown in Fig. 5(a). Each IFMD spectrum was measured by 
moving the slit (Slit 2) in front of the scintillation counter. 
Fig. 6 shows the IFMD spectra as a function of distance x 
parallel to the surface for values of D between 0 and —20 um 
(the minus sign indicates that the direction of the force is 
opposite to the force of gravity). When D = 0 um, neither a 
mirage diffraction beam nor IMFD are observed, so the 
crystal is not bent. When the absolute value of D increases, 
IFMD appear and the number of interference fringes 
increases. One, two, four and five peaks are observed for D = 
—5, —10, —15 and —20 urn, respectively. Similar variations 



have been reported by Jongsukswat et al. (2012a) in the 
Si(220) section topographies. The IFMD spectrum for D = 
0 um in Fig. 10 of their paper is quite similar to that for D = 
—10 um in the present Fig. 6. Note the different meaning of 
D = 0 urn in these two figures; D = 0 um in the former figure 
means that no force is applied at the free end, while D = 0 um 
in Fig. 6 of this paper means that the sample is flat. We can 
estimate the actual value of D in the former by assuming that 
the strain gradient is the same as that for D = — 10 um in the 
latter and the strain gradient is given as e' zz = 3vDl/L? under 
the cantilever load. In the experiment of Fig. 6, the original 
sample length is L = 40 mm, but it is likely that the effective 
length is shortened to 30 mm because of buckling. Then we 
obtain the effective displacement D elt = 26 um, while we 
obtain the value of Z> e(£ = 10 um if we use the original value of 
L = 40 mm. In either case, the effective displacement in the 
experiment by Jongsukswat et al. (2012a) is the same order of 
magnitude as that obtained in Fig. 4 (D = 28 um). 
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Figure 6 

Measured spectra of IFMD for / = 16 mm. D is varied as 0 um (filled 
circles), —5 |rm (open circles), —10 um (crosses), —15 um (triangles) and 
—20 um (squares). 
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4. Conclusions 

IFMD and IFMRB interference fringes were observed when 
no artificial force was applied to the sample crystal. The origin 
of the interference fringes is attributed to the strain induced 
by the uniformly distributed load due to gravity. According to 
rod theory, the strain gradient s' zz under a uniformly distrib- 
uted load is proportional to both I 2 and D and is given by 
equation (5), while it is proportional to both / and D under the 
cantilever load as given by e' zz — 3uD//L 3 . In the present 
experiment, the values of s' zz determined by measuring IFMD 
and IFMRB show I 2 dependence in accordance with rod 
theory under a uniformly distributed load, except for a region 
of large strain (and large I 2 ) near the fixed edge. IFMRB are 
observed even when I 2 = 0 mm 2 , at which point the strain 
gradient should be zero and no mirage diffraction should be 
observed. The observed IFMRB can be explained by the 
residual strain from manufacture. From the slope g t of the line 
I 2 versus e' zz in Fig. 4, we estimate the value of D to be 28 urn 
using equation (5). This value is of the same order of magni- 
tude as that estimated by IFMD under a cantilever load in the 
opposite direction to gravity. 

It should be noted that the IFMD and IFMRB interference 
fringes, accompanied by mirage diffraction, are very useful for 
measuring the strain distribution in a very weakly distorted 
crystal. The lower limit of detection of strain is of the order of 
1(T 7 . 

In addition, the measurement of strain distribution using 
IFMD and IFMRB has some advantages, i.e. the measurement 
can be done without directly touching the sample or 
destroying it, with high sensitivity, and at any depth from the 
top surface to the bottom surface. The present usage of IFMD 
and IFMRD can also be applied to a diffractometer with high 



energy resolution; this application will be investigated in the 
near future. 
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